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Rarita-Schwinger-Weyl semimetal in J e ff = 3/2 electron systems 
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We propose a relativistic J E // =3/2 semimetal with 4a ! 1 or 5 d 1 electrons on a cubic lattice when 
the strong spin-orbital coupling takes over the Hunds’ coupling. A relativistic spinor with spin 3/2 
is historically called Rarita-Schwinger spinor. In the massless case, the right- and left-handed chiral 
degrees of freedom of the Rarita-Schwinger spinors are independent. In the lattice model that we 
propose, the right- and left- handed gapless points in Brillouin zone are separated. We call this 
linearly dispersed semimetal Rarita-Schwinger-Weyl semimetal, similar to Weyl semimetal for spin 
1/2 systems. There is a network of gapless Fermi arcs in the surface Brillouin zone if ni ± Ji 2 ± «3 
is even for the normal vector (m, « 2 , n- 3 ) of the surface while the surface is insulator if ni + ri 2 + «3 
is odd. 


Introductions. In 1936, eight years after derived his 
famous relativistic electron equations of motion, Dirac 
generalized these equations to higher spin relativistic 
particles^. The first important example was the recov¬ 
ering of Maxwell equations. The next simplest particle 
except the electron and the photon is of spin-3/2 and 
obeys so-called Rarita-Schwinger (RS) equations^. This 
fermonic field later played an important role in super¬ 
gravity theory, known as the super-partner of the gravi¬ 
ton, the gravitinoi 

Electrons and photons are particles accompanying us 
on daily life while the RS particles are not found even 
in experiments of high energy physics or in cosmology 
observations. Recently, the interplaying between high 
energy physics and condensed matter physics supplies a 
new playground to the relativistic systems, e.g., Dirac 
semimetal in graphene^ and three dimensions^”—, topo¬ 
logical insulator^, supersymmetric systems^”—, and 
newly proposed—^ and discoveredLLSI Weyl semimetal. 
Can we expect a RS semimetal in this playground? 

The earlier theoretical prediction for Weyl semimetal 
was based on the J e ff = 1/2 states in iridatesi^ which 
belongs to a large class of 4 d and 5 d transition metal 
oxides with a strong spin-orbital coupling. When the d- 
orbitals are below 2/3 filling, the crystal field projects 
the electron configurations to the t- 2 9 orbitals. When the 
strong spin-orbital coupling dominates, the electrons fill 
the J e ff = 1/2 or J e ff =3/2 depending on the filling fac¬ 
tor. For iridates, the 5 d 5 electrons occupy all J e // = 3/2 
states and half fill J e // = 1/2 orbitals^. For d 1,2 elec¬ 
trons, they quarterly or half occupy the J e // = 3/2 
states. The representatives of these J e // = 3/2 materials 
are the ordered double perovskites with the chemical for¬ 
mula A^B'BOq where B' ions are commonly 4d 1,2 and 
5d 1,2 transition metals’, e.g. Mo +5 ,Re +6 ,Os +7 for d 1 and 
Re +5 ,Os +6 for d 2 . Several exotic magnetic phases and a 
quadrupolar phase were presented in these strongly cor¬ 
related materials^”—. 

Due to the particle-hole symmetry, the relativistic 
physics usually emerges in a half-filled lattice model. For 


the massive RS theory, this implies a 4eP or 5 d 2 system 
is required. However, the above strongly correlated sys¬ 
tems are not likely because the first Hund’s rule due to 
the interaction makes two electrons first form a total spin 
S = 1 state which then couples to 1 = 1 state. The final 
effective state is projected to J e // = 22i. On the field 
theoretical side, there are also many complexes for the 
massive RS field itself. We thus restrict our attention to 
the massless RS field. 

The massless relativistic spinor systems with any spin 
possess only two physical degrees of freedom, the high¬ 
est and lowest helical states^. Maxwell field with only 
transverse components is the well-known example. The 
solutions of the massless RS equations possess the helicity 
±3/2 associated with the right- and left-handed chirality. 
The helicity ±1/2 states are not zero mass states. There¬ 
fore, the fourfold-degeneracy of the J e // = 3/2 bands in 
the transition metal compound where the massless RS 
field emerges needs to be lift: Two helicity ±1/2 bands 
are gapped from the helicity ±3/2 bands. The relativistic 
theory emerges from the half-filling of two lower bands, 
i.e., the quarter filling J e ff = 3/2. We then see that ap¬ 
propriate systems are 4c! 1 or 5 d 1 transition metal com¬ 
pounds. They are semi metallic, similar to the Dirac and 
Weyl semimetals. We call these semi-metallic states the 
RS or RS-Weyl semimetal. We will present a simplest 
lattice model where the RS field emerges while searching 
for realistic materials is not our goal at this stage. 

On the other hand, we will not concern the Lagrangian 
presented in Rarita and Schwinger’s original paper in 
1941—. A modern version for the gravitino will be used 
because of its simplicity and gauge invariance^. The lat¬ 
ter is also that high energy theorists are studying and 
well-behaved in an external fields. 

To our knowledge, the RS physics was not concerned in 
the condensed matter physics context. We therefore will 
briefly introduce the RS theory and make the conven¬ 
tion to be familiar with the readers in condensed matter 
physics area. We then construct a tight-binding model 
for d xy , d yz and d zx electrons on a cubic lattice where 
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d xy refers to the d-e lectron annihilation operator with 
the a;j/-orbital component, etc. The long wave length ef¬ 
fective theory shows that this is a RS-Weyl semimetal 
with eight chiral Weyl points. The surface states are de¬ 
pendent on what surface of the system is. If the normal 
vector ( 711 , 712 , 713 ) of the surface Brillouin zone is even, 
i.e., 7 ii + 7 i 2 + 713 is even, there is a network of the gap¬ 
less Fermi arcs. For a surface Brillouin zone with an odd 
normal vector, the surface is insulating. The example of 
the former is the surface [ 110 ] while [ 001 ] and [ 111 ] are 
the representatives of the latter. 

Rarita-Schwinger theory. The RS equations for a sixteen 
component vector-spinor field t/> mq in 3+1 dimensions are 
given by 


~ m )^ = 0, (1) 

X = 7^ = 0, (2) 


where the convention we use are: /i = 0,1,2,3 denote 
the time-space indices with flat metric 77 °° = 7700 = 1 and 
1 in = — 77 " = —1; a = sa for s = R, L and a =t,i are 
chiral and spin indices, respectively. The Weyl’s gamma 
matrices are chosen: 
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In Eq. ©, the four vector indices are contracted over 
so that x is a pure Dirac spinor, x = 0 projects out the 
spin- 1/2 sector and leaves only the degrees of freedom 
of the spin-3/2 sector. A non-zero mass m brings many 
problems to the RS system. It was known that if m ^ 0, 
there will be fermonic modes with superluminal velocities 
if the RS field couples to the external electromagnetic 
field in a minimal wa y 27 ' 28 . However, this does not cause 
a trouble for the condensed matter system because it is 
allowed that the velocity of a collective mode exceeds the 
’speed of light”, the Fermi velocity. 

The bad thing for us is that the four vector leads to 
a superfluous component 7/0 if we identify i = 1,2,3 
to the orbital degrees of freedom xy,yz, zx of the d- 
electron. One may use Eq. © to eliminates 7 / 0 . The 
canonical quantization procedure is very inconvenient be¬ 
cause the Hamiltonian after eliminating 7/0 becomes very 
complicated^, although a quantized free RS field is well 
defined^. Thus, it is difficult to realize a massive RS field 
with a simple model. We will study a massless RS field 
in this work. We take a simple form which used in super¬ 
gravity for the gravitino RS field. In a flat space-time, 
the RS Lagrangian reads^ 

£ = (3) 

where 7 /+ = ippl 0 and 7 y/iA = i s the total an- 

tisymmetricized product of three gamma matrices. It 
is obvious that the Lagrangian is gauge invariant under 
Vv — > ipii + d^e for an arbitrary spinor e. The equations 
of motion can be written in many equivalent ways^, e.g, 

rd^ v - d vX = 0. (4) 


in terms of Lagrangian ©. Imposing the gauge fixing 
X = 0, Eq.© recovers the massless RS equations. It is 
also easy to have a supplementary condition = 0 

according to the RS equations. 

7 / 0=0 gauge and two component chiral spinors. Since 
7 is four component vector, it is still difficult to con¬ 
nect with a condescend system. The gauge invariance of 
Lagrangian m enables us to take t/ 0 = 0 gauge as tak¬ 
ing Aq = 0 gauge for the electromagnetic field. In this 
gauge, the supplementary condition are diip 1 = 0, which 
also projects the system to S = 3/2. 

To more explicitly link to a condensed matter model, 
we use two component forms of the chiral RS spinor , 
i.e., 7 l> ia = (ciLa ) cm a ) T . Making Fourier transforma¬ 
tions c iR ( L ) a = J2 P c R(L)ape lPlJ - xl * with py = (w, p), the 
Lagrangian for 7/0 = 0 is given by 


£ 


uj.p 


^ ^ v C isp C sp 
s=L,R 

£ ijk 4L P Pj c kLp - 3k c\ RpPj c kRp . 


(5) 


This Lagrangian is separated into the right- and left- 
handed chiral sectors as that for the massless spin- 
1/2 Dirac spinor. The supplementary conditions are 
P l Ci sp = 0 . is a residual gauge invariance as Ci S (x) —> 
Cis(x) + 9ie(r) if the spinor e(r) is independent of time. 
To fix the gauge invariance, we impose a l Ci sp = 0. The 
Hamiltonian then can be read out 


n p = v 3/2 ni n v 3/2 , ( 6 ) 

K n = - ieijk 4LpPj c kLp + ie ijk cl RpPj c kRp , (7) 

where V 3 / 2 denotes the projection to the spin-3/2 sector. 
The matrices of the left- and right handed unprojected 
Hamiltonians 7are of a simple form 


^+L;6x6 = ±P-T®/ 2 (8) 

where (T i y k = ie l i k are three generators of SO(3) group 
and I 2 is 2 x 2 unitary matrix for spin index. Diagonaliz¬ 
ing 7 i. can be done by diagonalizing 77/ ra and making an 
projection V 3 / 2 . It is easy to see the eigen values of 7 V^ R 
are w = 0 and w = ±|p[. The supplementary condition 
suppresses the eigen state corresponding to w = 0. The 
solutions of u> = ±|p| are not independent: The positive 
(negative) solution for the right-(left-)handed RS field is 
particle while the other is hole,or antiparticle. It is easy 
to check the right-(left-)handed solution is of the helicity 
3/2 (-3/2)25^. Therefore, projecting to J e // = 3/2, the 
genuine physical degrees of freedom for the ground states 
are 1+1. The electron operators can be expanded as 


Ci R (L)a — y^jc ±3 / 2j) U R ( L ) ia e zpx + d± 3 / 2 ,p' |/ fl(L)i (T e lpx ) 
p 


with a l U R (L)i = <J l V R (L)j = 0. The solution vector- 
spinors U and V are normalized by U)] o 3 U lR = 
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V l ^a j ViR = —pi/u and U R ( L ^iR{L) = 0. U and V are 
not independent and they are related by the charge con¬ 
jugation. 

With the RS field operators c± 3 / 2 , P and d ± 3 ^. 2 p in the 
diagonalized basis, 

CiR(L)a = y^Xc±3/2,p^R(L)iae lpX + <^± 3 / 2 ^R( L )ia e ~ lPX ) 

P 


with a l Uji^L)i = o l Vfi(L)i = 0 . The diagonalized Hamil¬ 
tonian then reads 

H = IPl(4p 5 Ap + 4pdAp)- (9) 

p,A=±3/2 


Defining a unit vector n = p/|p| according to (0 and 
after the projection, a topological winding number is ex¬ 
plored for a chiral RS field 


vr,l 



l2 dn dn 
d pn ■ -— x -— 

&Px ^Py 


( 10 ) 


where S 2 is a two-sphere which wraps the Weyl point in 
the wave vector space. 

Lattice model and RS-Weyl semimetal . We now turn 
to the realization of RS field in condensed matter sys¬ 
tem. As we mentioned, we assume the 4 d 1 or 5c? 1 
ions are on a simple cubic lattice for transition metal 
(e.g., Mo +5 ,Re +6 , Os +7 ) compounds. The crystal field 
projects the effective orbital angular momentum to 1 = 1. 
We consider the strong spin-orbital coupling takes over 
the Hunds’ coupling, the d 1 electrons fill one of four 
Jeff = 3/2 states, i.e., the quarter-filling of J e ff = 3/2. 
The Hamiltonian we are studying is as follows 

a 

± d\z,adxy,a±S J + h.C. - A ^ l a • S a 4-(11) 


where S xyz are the lattice vectors in the positive di¬ 
rections. The hopping term is between the different 
orbitals in different direction. Furthermore, the hop¬ 
ping in the negative Si direction carries a phase 7r while 
carrying no phase in the positive direction (See Fig IT]). 
The third term is the spin-orbital coupling projected to 
Jeff = 3/2^2. The ■ ■ • refers to other hopping terms like 

td 'Ei=xy,yz,zx YX(ab) d\adib- TVe assume td is small and 
neglect this term for the moment. 

It is easy to check that the Hamiltonian is not invari¬ 
ant under space inversion P: Under di a —> —di t - a , t 
term changes sign. However, it is time reversal T in¬ 
variant: Under d iR ( L)<J , a -)■ ad iR ^_ a ^ H H. Such 
properties under P and T resemble the recent observed 
Weyl semimetal in the non-centrosymmetric transition 
metal monophosphidesi^— . The role of the strong spin- 
orbital term is separating the J e // = 1/2 orbitals from 
J e f f = 3/2 ones with an order of leV gap. Due to the mi¬ 
nus sign before the spin-orbital coupling, the J e // = 3/2 


states have a lower energy. The effective Hamiltonian in 
the wave vector space is given by 

H = V 3/2 H un V 3/2 , (12) 

H un = -2it^2e ijk smkjd^dku 

k 

where the lattice spacing is taken to be one. Therefore, 
k = 0 and k = it are Weyl points for the left- and right 
handed chiral RS spinors by denoting dqo+< 5 k = CiL, P 
and = Ci Rp . Hamiltonian (fl2l) exactly recovers 

the RS Hamiltonian 0 with vf = 2t = 1. The helicity 
±1/2 states are projected over from the massless sector. 
Thus, Hamiltonian CD describes Rarita-Schwinger- Weyl 
semimetal. 


0 Jt 


+ < - 
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FIG. 1: (color online) The hopping of electrons with different 
orbital components. There is a 7r phase when electron hops 
in the negative direction. 


In fact, there are other three pairs of Weyl points 
which are located {(0, n, n), (n, 0,0)}, {(7r, 0,7r), (0, n , 0)} 
and {(7r, 0,0), (0,0,7r)} beside of the pair {0, fr}. In these 
points, there are singularities of the vector field 


sink 


(13) 


Corresponding to these singularities, the winding num¬ 
bers defined as Eq. m are (1,-1) for each pair, respec¬ 
tively (See Fig. [2]). These Weyl points can be thought of 
as the ” monopoles” in Brillouin zone whose configuration 
is scratched in Fig. [2] 

Surface states and Fermi arc network. One of the char¬ 
acteristics of the Weyl semimetal for spin-1/2 system is 
the existence of Fermi ares in surface states^, which also 
has been observed experimentally. We now study the 
surface state of the RS-Weyl semimetal. We consider a 
surface Brillouin zone with the normal vector along the 
(ni,n 2 , n 3 ) direction. The bulk Weyl points are projected 
onto this surface. A set of the Weyl points may have the 
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FIG. 2: The singularity (monopole) configurations and dis¬ 
tributions, as well as the flows of the flux in Brillouin zone. 
The eight monopoles are located at (0, 0, 0), (7r, 0, 0), (0, n, 0), 
(0,0,7r), (-7T, 7r,0), (tt, 0,7r), (0,7r, 7r) and (tt, with the 
winding numbers (1, —1, —1, —1,1,1,1, —1). 


same image. The charge of the image is defined by the 
sum of the winding numbers of these Weyl points in this 
set divided by the number of the points. It is easy to see 
the charges of all images are zero if ni + ri 2 + n 3 is odd 
while the charges of the images are ±1 if n\ + 712 + 713 
is even. For a given surface, these images form a lat¬ 
tice. For a lattice where all sites are zero charged, there 
is no reason to see the Fermi arc connecting any pair 
of images while the Fermi arcs may appear between the 
positive and negative charged images in a charged site 
lattice. In Fig. EKa)), we sketch the image network for n 
even surface, the [110] surface. The black and while spots 
correspond to the images with charge ±1. We numeri¬ 
cally calculate the surface states of the [ 001 ], [ 111 ] and 
[110] surfaces. We see that there is no gapless state for 
[ 001 ] and [ 111 ] while the dispersions of the surface states 
for the [110] are shown in Fig. Efb). For the latter, as 
expected, the gapless Fermi arcs form a rectangular net¬ 
work and the images are the lattice sites of the network. 

We now analytically verify the insulating nature of the 
odd surface. For example, we can show the insulating 
state in the [ 001 ] surface which consists of two most-up 
coupled [001] layers, (x,y,L/ 2) and (x, y, L/2 — S z ). The 
effective surface Hamiltonian is given by 


#001 = # 3 / 2 # s “ n # 3/2 

#001 = ~t £ (±4„ ;Oia d ax;o±a# , 0 


(14) 


a x ,a v :cx.—ua 


d , y Z ;a I ad'Zx;a±5 x ,a) 

^-L (—dyZ'a udxy,a,d 

OX lOy 

+ dyz t a,ddxy,a, u ) + h.c. 
td d iaa diba, 

i,ct,(ab) 


(15) 



FIG. 3: (color online) The Fermi arc network of the surface 
states in the [110]-surface Brillouin zone. Upper: Images of 
the Weyl points with charges ±1 (black and white). Lower: 
The dispersions and the network of gapless Fermi arcs for 
t = 1 . 


where t± is coupling between two layers and determined 
by the boundary condition; a = u,d are layer indices. 
We include the diagonal hopping term with the hopping 
amplitude td which may affect the surface state. One 
eigen value of #qo\ is uj = —td(cosk x + cos k y ) which is 
projected out after T 3/2 projection. The real spectra of 
#001 are given by 

uj = — td(cosk x + cosky) (16) 

± 2 y£ + t 2 sin 2 k x + t 2 sin 2 k y . 

We assume td < t. When t± > td, the surface state is 
insulating which is indeed that as we expected. When 
t _l = td, the surface becomes a semimetal with quadratic 
dispersion at (0,0) and (n, n). If t± < td, two Fermi 
pockets appear at (0,0) and ( 7 r, 7 r). Surrounding (0,0), 
it is a left-handed helical particle Fermi pokiest while 
another is hole’s with a right helicity surrounding ( 7 r, tt) 
(See Fig. H) 

We can also show that there is the Fermi arc network 
for an even surface by a similar analysis. 

Conclusions and perspectives. The spin -3/2 relativistic 
spinor state, the RS semimetal, was proposed in tran¬ 
sition metal compounds with J e ff =3/2 4 d 1 and hd 1 
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FIG. 4: (color online) The Fermi pockets of the surface states 
in the [001]-surface Brillouin zone for td > t±. The blue 
pocket is particle’s and the red one is hole’s. 

electron configurations. The fruitful surface state prop¬ 
erties were discussed and the network of the gapless Fermi 
arcs was predicted. Our results point a new direction for 
the researches of topological states of matter after the 


topological states of spin-1/2 Dirac, Weyl and Majorana 
fermions are exhaustively studied. 

Unlike the spin-1/2 spinor whose properties are very 
well understood in field theory, many properties of spin- 
3/2 relativistic field are not transparent. We were not 
aware of the topological classification of spin-3/2 spinor 
although the index theorem for the RS operator was 
known. The massive RS theory may also be correspond¬ 
ing to the topological insulator while the lack of gauge 
invariance leads to the Hamiltonian of the system be¬ 
comes very complicated. The existence of the fermionic 
modes propagating faster than ’’the speed of light” in 
the minimal coupling to the external electromagnet field 
is expected to be experimentally examined. The topo¬ 
logical superconductivity of spin-3/2 RS spinors were not 
studied. We also anticipate the realization of the RS field 
in cold atom systems. 
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